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Instability in the Gel'fand inverse problem at 

high energies 



M.I. Isaev 



f^i . Abstract 

We give an instability estimate for the Gel'fand inverse boundary value 
£^ ■ problem at high energies. Our instability estimate shows an optimality of 

^ ' several important preceeding stability results on inverse problems of such 

a type. 



1 Introduction 

In this paper we continue studies on the Gel'fand inverse boundary value prob- 
lem for the Schrodinger equation 

— Atp + v(x)ijj — Etp, xeD, (1.1) 

where 



D is an open bounded domain in Mr, d > 2, 
with dD e C 2 , 



(1.2) 



<X>. veh°°(D). (1.3) 



As boundary data we consider the map <3E> = <&(E) such that 



Hm^lav) = ^\8D (1-4) 

o 

for all sufficiently regular solutions ip of (jl.ip in D = D U dD, where v is the 
outward normal to dD. Here we assume also that 

E is not a Dirichlct eigenvalue for operator —A + v in D. (1-5) 

X' 

The map $ = 4>(E) is known as the Dirichlet-to- Neumann map. 

We consider the following inverse boundary value problem for equation (jl.ip : 

Problem 1.1. Given $ for some fixed E, find v. 

This problem is known as the Gel'fand inverse boundary value problem for 
the Schrodinger equation at fixed energy (see [7], [!!])■ At zero energy this 
problem can be considered also as a generalization of the Calderon problem of 
the electrical impedance tomography (see [S], Q3Ij). Problem 1.1 can be also 
considered as an example of ill-posed problem: see [2], [3] for an introduction 
to this theory. 

There is a wide literature on the Gel'fand inverse problem at fixed energy. 
In a similar way with many other inverse problems, Problem 1.1 includes, in 
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particular, the following questions: (a) uniqueness, (b) reconstruction, (c) sta- 
bility. 

Global uniqueness results and global reconstruction methods for Problem 
1.1 were obtained for the first time in |19] in dimension d > 3 and in [3J in 
dimension d = 2. 

Global logarithmic stability estimates for Problem 1.1 were obtained for the 
first time in [JJ in dimension d > 3 and in [25 in dimension d = 2. A principal 
improvement of the result of [JJ was obtained recently in [23] (for the zero energy 
case): stability of [24] optimally increases with increasing regularity of v. 

Note that for the Calderon problem (of the electrical impedance tomography) 
in its initial formulation the global uniqueness was firstly proved in [30) for d > 3 
and in [T? for d = 2. Global logarithmic stability estimates for this problem 
were obtained for the first time in [JJ for d > 3 and [T5J for d = 2. Principal 
increasing of global stability of [JJ , [T5J for the regular coefficient case was found 
in [21] for d > 3 and [25] for d = 2. In addition, for the case of piecewise real 
analytic conductivity the first uniqueness results for the Calderon problem in 
dimension d > 2 were given in [13] . Lipschitz stability estimate for the case 
of piecewise constant conductivity was obtained in [2J (see [27J for additional 
studies in this direction). 

The optimality of the logarithmic stability results of [JJ , [IS] with their princi- 
pal effectivizations of [24], [28] (up to the value of the exponent) follows from [16] . 
An extention of the instability estimates of [TB] to the case of the non-zero en- 
ergy as well as to the case of Dirichlet-to-Neumann map given on the energy 
intervals was obtained in [5J. 

On the other hand, it was found in [50], [2JJ (see also [33], [35]) that for 
inverse problems for the Schrodinger equation at fixed energy E in dimension 
d > 2 (like Problem 1.1) there is a Holder stability modulo an error term 
rapidly decaying as E — > +oo (at least for the regular coefficient case). In 
addition, for Problem 1.1 for d = 3, global energy dependent stability estimates 
changing from logarithmic type to Holder type for high energies were obtained 
in [13], [UJ. However, there is no efficient stability increasing with respect 
to increasing coefficient regularity in the results of [12]. An additional study, 
motivated by [T2J, [23], was given in [T5] . 

The following stability estimate for Problem 1 . 1 was recently proved in [JJJ : 

Theorem 1.1 (of 11). Let D satisfy 1(13)1 . where d>3. Let Vj G W m ' x {D), 
m > d, supp Vj C D and \\Vj\\w m -i-(D) < N for some N > 0, j — 1,2, (where 
W m ' p denotes the Sobolev space of m -times smooth functions in ]L P ). Let v±, «2 
satisfy (|1 .5)> for some fixed E > 0. Let <f>i(E) and $2(-E) denote the DtN maps 
for vi and vi, respectively. Let s\ = (m — d)/d. Then, for any r G (0, 1) and 
any a, j3 G [0, si], a + j3 — s±, 

\\v 2 - ui||ioo (D) < A(l + VE)6 T + B(l + VE)~ a (In (3 + (T 1 ))^ , (1.6) 

where S — \\<£>2{E) — < 1 > i(£')||l° c (9D)-!.l oo (9D) an d constants A, B > depend 
only on N , D, m, r. 

In particular cases, Holder-logarithmic stability estimate (11.61) becomes co- 
herent (although less strong) with respect to results of [21], [23], [24]. In this 
connection we refer to [TT for more detailed infromation. Concerning two- 
dimensional analogs of results of Theorem 11.11 see [2D] , [2SJ , [25J , [25] . 
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In a similar way with results of [DJ, [TDJ, estimate (ll.6[) can be extended to 
the case when we do not assume that condition (|1.5[) is fulfiled and consider an 
appropriate impedance boundary map (or Robin-to- Robin map) instead of the 
Dirichlet-to-Neumann map. 

In the present work we prove optimality of estimate (|1.6|) (up to the values 
of the exponents a, /3) in dimension d > 2. Our related instability results for 
Problem 1.1 are presented in Section 2, see Theorem 12.11 and Proposition l2.il 
Their proofs are given in Section 4 and are based on properties of solutions of 
the Schrodinger equation in the unit ball given in Section 3. 

2 Main results 

In what follows we fix D — B d (0, 1), where 

B d {x°,p) = {x eR d :\\x- x a \\ Ed < p}, x eR d ,p>0. (2.1) 

Let 

|I_F|| denote the norm of an operator 

11 11 (2.2 

We recall that if v\, v 2 are potentials satisfying (|1.3I) . (| 1 . 5|) for some fixed E, 
then 

§ 2 {E) - &\{E) is a compact operator in L°°(&D), (2.3) 

where $i, <3E>2 are the DtN maps for Vi, v 2 , respectively, see [TD], [22]. 
Our main result is the following theorem: 

Theorem 2.1. Let D = B d (0, 1), where d > 2. Then for any fixed constants 
A, B, k,t,s > 0, m > d and s 2 > m there are some energy level E > and 
some potential v € C m (D) such that condition (|1.5[) holds for potentials v and 
vq = 0, simultaneously, supp v C D, ||w||l°°(_D) < £> || w ||c m (-D) — ^> w ^ ere 
Ci = Ci(d,m) > 0, 6ut 

||«-«b|U- (1J ) > A(l + VE) K S T + B(1 + VE) 2 ^ S - S ^ + S- 1 )y s (2.4) 

/or any s € [0, S2], where $0 are DtN map for v and vq, respectively, and 
5 = \\&(E) - $o(-E)|| is defined according to (j2~2j) . 

Theorem 12.11 shows . in particular, the optimality (at least for potentials in 
the neighborhood of zero) of estimate (|1.6[) (up to the values of the exponents a, 
/?). As a corollary of Theorem 12. 1[ one can obtain an optimality of the stability 
results of [2D], [2T], [23], [25]. 

In the present work Theorem 12.11 is proved by explicit instability example 
with complex potentials. Examples of this type were considered for the first 
time in [16] for showing the exponential instability in Problem 1.1 in the zero 
energy case. An extention to the case of the non-zero energy as well as to the 
case of Dirichlet-to-Neumann map given on the energy intervals was obtained 
in[S]. 

Let us consider the cylindrical variables: 

(ri, 9, x') £ R+ x R/2ttZ x R d ~ 2 , 

ricos0 = xi, n sin 9 = x 2 , (2-5) 
x' = (x 3 , ■ ■ .,x d ). 
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Take </> G C°°(R 2 ) with support in B 2 (0, l/3)n{x x > 1/4} and with ||<£|| L =o = 1. 
For integers m,n > 0, dehne the complex potential 

v nm = n- m e in U(ri,\x'\). (2.6) 

We recall that 

||w„ m ||L=° = vT m , ||u m „||c™ < Ci, (2.7) 

where Ci = G\(d,m) > 0. Note that C% is the same as in Theorem 12.11 
Estimates (|2.7p were given in [TB] (see Theorem 2 of |16j). 

To prove Theorem 12. II we use, in partucular, the following proposition: 

Proposition 2.1. Let D = B d (0, 1), wh ere d > 2. Lei condition (|1.5[) /io/d wi</i 
i> = w„ m (o/ (|2.6[) J and u = i>o = for some E > and some integers m > 0, 
n > 20(1 + \/S) 2 . T/ien. /or any a > 0, 

\\$ nm (E) - M e )\\h-"(s^)^(s^) <C 2 (1 + Q + EQ)2- n '\ (2.8) 

where $ nm , $o i/ie DtN map for v nrn andvo, respectively, C2 — C2(d,a) > 0, 

Q = II (- A + v o - ^) _1 Hl2(d)^l2(d) + II (-A + v nm - Ey 1 \\ h 2(D)^(D), (2.9) 

where (—A + vq — E)~ l , (—A + v nm — E)^ 1 are considered with the Dirichlet 
boundary condition in D and H ±rT = W ±a ' 2 denote the standart Sobolev spaces. 

Analogs of estimate f|2.8[) (but without dependence of the energy) were given 
in Theorem 2 of [TB] for the zero energy case and in Theorem 2.4 of [S] for the 
case of the non-zero energy and the case of the energy intervals. 

We obtain Theorem 12.11 combining known results on the spectrum of the 
Laplace operator in the unit ball (see formula (|4.9[) below), Proposition 2.1, 
estimates (|2.7|l and the fact that 

ll^llL°°(s d - 1 )^L~(s t! - 1 ) < c(d, a)\\F\\ H -v( S d-i)^ H v( S d-i) (2.10) 

for sufficiently large a. The detailed proof of Theorem 12.11 and the proof of 
Proposition are given in Section 4. These proofs use, in particular, results, 
presented in Section 3. 

Remark 2.1. In a similar way with [TB], [5], using a ball packing and covering 
by ball arguments (see also [B]), the instability result of Theorem 12.11 can be 
extended to the case when only real-valued potentials are considered and in the 
neighborhood of any potential (not only vq = 0). 

3 Some properties of solutions of the Schrodinger 
equation in the unit ball 

In this section we continue assume that D = B d (0, 1), where d > 2. We fix an 
orthonormal basis in JL 2 ^- 1 ) = L 2 (<9£>) 

{fj P -j> 0, 1 <p<Pj], , 3 
fj p is a spherical harmonic of degree j , 
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where pj is the dimension of the space of spherical harmonics of order j, 



where 



and 



n\ n(n — 1) • • • (n — k + 1) , 
] —-t for n > (3.3) 



k k\ 



for n < 0. (3.4) 



The precise choice of fj P is irrelevant for our purposes. Besides orthonormality, 
we only need fj P to be the restriction of a homogeneous harmonic polynomial of 
degree j to the sphere and so \x\i fj P {x/\x\) is harmonic. We use also the polar 
coordinates (r,w) e R + x § d_1 , with x = ruj € R d . 

Lemma 3.1. Let D = B d (0, 1), where d > 2. Let potential v satisfy (|1.3[) and 
(|1.5[) for some fixed E. Let ||i>||jl«>(.d) < N, for some N > 0. Then for any 
solution ijj € C(D U dD) of equation the following inequality holds: 

uw^-(d) < (i + (jv+ i^i)ii(-A + v - sr 1 !!^^^)) \\f\y {d D), (3.5) 

where f = tp\dD> ( — A + w — -E) -1 * s considered with the Dirichlet boundary 
condition in D. 

Proof of Lemma \3.1\ We expand the function / in the basis 

f = J2 c jpfj P - ( 3 - 6 ) 



We have that 



Let 



-Em 2 - ( 3 - 7 ) 



Il 2 (o_d) — 



Note that 



M x ) = E c Jp r " '/j>M- ( 3 - 8 ) 
IIV'oIIl^d) = E I c ip 1 2 1 1 ^ /jp (^) 1 1 2.2 (j? ) = 

,i (3.9) 

= Em 2 / ^ +!i ^<Ei^i 2 

Jo 

Using (jl.lj) and the fact that V'o is harmonic, we get that 

(-A + v-E){i}-^ ) = (E-v)i> . (3.10) 
Since ip\dD = ipo\dD = f, using (|3 . 10[) . we find that 

U - MupD) < (N +\E\)\\(-A + v- Ey^^^^W^oW^^). (3.11) 
Combining $£7$, ([3T9]) . (|3~TT|) . we obtain ([33]) . ■ 
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Let < •, • > denote the scalar product in the Hilbert space h 2 (dD): 



<f,9>= J f(x)g{x)dx, (3.12) 

dD 



where f,g G l?(dD). 



Lemma 3.2. Let D = B d (0,l), where d > 2. Let potentials v\, v 2 satisfy 
(|1.3|) and (| 1 - 5[) for some fixed E. Let v\, vi be supported in B d (0,l/3) and 
\\ v i\\h°°(D) < N, i = 1,2, for some N > 0. Then for any ji,j2 £ HU {0}, 

1 < Pi < Pn, 1 < Pi < Pj 2 and j max = max{ji,j 2 } > 10(1 + VT^I) 2 the 
following inequality holds: 

(f jlPl , ($i(S) - < C(d)(l + (JV+|JS|)0)2-^~, (3.13) 

w/iere 

g = IK-A + W! - J B)- 1 || L 2 (I?) ^ L 2 (D) + ||(-A + w 2 -S)- 1 || L 2 (D) ^ L 2 (I 3 ) , (3.14) 

$i, $2 are i/ie Z)iA^ map for v\ and v^, respectively, and (—A + v% — E)^ 1 , 
(— A + i?2 — -E) -1 are considered with the Dirichlet boundary condition in D. 

Analogs of estimate (13. 13)) (but without dependence of the energy) were 
given in Lemma 1 of |16j for the zero energy case and in Lemma 3.4 of [5] for 
the case of the non-zero energy and the case of the energy intervals. 

We prove Lemma \'S. 2 1 for E ^ in Section 5, using expression of solutions 
of equation —Atp = Eijj in B d (0, 1) \ B d (0, 1/3) in terms of the Bessel functions 
J a and Y a with integer or half-integer order a. 



4 Proofs of Proposition 2.1 and Theorem 12.1 



We continue to assume that D = B d (0, 1), where d > 2 and to use the orthonor- 
mal basis {f jp : j e N U {0}, 1 < p < pj} in L 2 ^" 1 ) = L 2 (&D). The Sobolev 
spaces H (7 (E> d ~ 1 ) can be defined by 



^^ c jpfjp '■ c jp/jp| < 

,p hv 

1 1 2 



(4.1) 



2c7| |2 



see, for example, [IB]. 

Consider an operator A : i? _<T (S <J_1 ) — > H a (§ d ^ 1 ). We denote its matrix 
elements in the basis {fjp} by 

a jipij2P2 fjipii Afj 2 p 2 > • (4-2) 

We identify in the sequel an operator A with its matrix {aj lPl j 2P2 }. In this 
section we always assume that ji, j 2 € N U {0}, 1 < pi < , 1 < p 2 < p j2 . 
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We recall that (see formula (12) of |16| ) 

PIIh—OS*- 1 )-*-*^^- 1 ) < 4 sup (1 + max{ji, j 2 }) 2<T+ ' i |aj 1 p 1 j 2 p 2 |- (4.3) 

Proof of Proposition 2.1. In a similar way with the proof of Theorem 2 of |16j 
we obtain that 

< f jlPl , (knn(E) - * {E)) f nP2 >= (4.4) 

for jmax — rnaxjjii J2} < [^^T^l (^ ne om y difference is that instead of the 
operator —A we consider the operator —A — E), where [•] denotes the integer 
part of a number. Note that 

71-1 



t + 1 > n/2 > 10(1 + V£) 2 , Km||L-(D) < I- (4-5) 
Combining (|Q|l . (|4~4)) . (|4~5l) and Lemma IO we get that 

||*mn(-E)-$o(-E)||//-''(S'i-i)^H-"(Sd-i) < 

< 4C(«Q ( 1 + (1 + £)q) sup (1 + i maa ) 2CT+d 2-^« < (4 6 ) 

V ' j ma *>n/2 V ' 

<c 2 (d,a)(i + g + sg)2-"/ 4 , 

where 

Q = \\(-A + v - £ , ) _1 ||l2(d)^l2(d) + ||(-A + w nm --B)" 1 || L 2 (I) )^ ]L 2( £ ,). (4.7) 



Let N(p) denote the counting function of the Laplace operator in D 

N(p) = |{A < p 2 : A is a Dirichlet eigenvalue of -A in D}\, (4.8) 

where | • | is the cardinality of the corresponding set. We recall that according 
to the Weyl formula (of |31j): 

N(p) < Cl (d)p d . (4.9) 

Lemma 4.1. Let D = B d (0, 1), where d > 1. Then for any p > 1 there is some 
E = E(p) G (p 2 ,2p 2 ) such that the interval 

(E(p)-c 2 p 2 - d ,E(p)+c 2 p 2 - d ) (4.10) 

does not contain Dirichlet eigenvalues of —A in D, where c 2 = 02(d) > 0. 

Proof of Lemma 4.1. We put c 2 = 2 d ~ 1 /(c 1 (d) + 1). Then we can select k 
disjoint intervals of the length 2c 2 p 2 ~ d in the interval (p 2 , 2p 2 ), where 



k = 



P 2 



2c 2 p 



2-d 



= [(c 1 (d) + l)p d ]>N(p). (4.11) 



Thus, we have that at least one of these intervals does not contain Dirichlet 
eigenvalues of — A in D — B d (0, 1). ■ 
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Proof of Theorem \2.1\ Let E = E(p) be the number of Lemma 4.1 for some 
p > 1. Using (I4.10[) . we find that the distance from E to the Dirichlet spectrum 
of the operator —A in D is not less than c 2 p 2 ~ d . Using also that E € (p 2 ,2p 2 ), 
we get that 

|| (-A - £0 _1 ||w(u)->^(D) < < E (d - 2)/2 /c 2 , (4.12) 

where (— A — E) is considered with the Dirichlet boundary condition in D. 
Let _ 

n = [20(1 + VE) 2 } + 1. (4.13) 

Using (|2.7[) and (|4.10p . we find that the distance from E to the Dirichlet spec- 
trum of the operator —A + v nm in D is not less than C2p 2 ~ d — n~ m , where v nm 
is defined according to (|2.6p . Since m > d and E e (p 2 ,2p 2 ) 7 using (|4. 13[) . we 
get that 

||(-A + v nm - 1 ||l 2 (d)— >l 2 (_d) < c 3 E (cl - 2 ^ 2 , 

E = E(p), p> Pl (d,m)>l, (4.14) 
C3 = c 3 («f,m) > 0, 

where (—A + t> nm — E)^ 1 is considered with the Dirichlet boundary condition 
in D. 

Combining Proposition H?T] and estimates (|2.10[) . (|4.12l) . (|4.14|) . we find that 

5 = \\$nm{E) - io(£)|k~(S<*-i)^L»(S^- 1 ) < Ci E d ' 2 2- n '\ 

E = E(p), p> Pl {d,m)>l, 

(4.15) 

n = [20(1 + VE) 2 } + 1 
C4 = C4(d, m) > 0. 

Since s 2 > m, taking p big enough and using (|4.15p . we obtain the following 
inequalities: 

n- ,n < e, (4.16) 
A(l + VE) K 5 T < \n~ m , (4.17) 

B(l + VE) 2 ^ (In (3 + S- 1 ))- 8 < \n~ m , 

< s < s 2 , 

where 

E = E(p), n = [20(1 + VS) 2 ] + 1. (4.19) 
Combining (J2U), (JH7J) , (l4TTo) - (j4~T^l) . we get that 

A(l + VE) K S T + B(l + VE) 2 ^-^ (In (3 + <T 1 ))~ S < 

||«nm||]L«>(D) = < £, 

|Pnm||c"'(-D) < Cl> 
supp v m „ C D. 



(4.20) 
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5 Proof of Lemma 13.21 



To prove Lemma 13.21 we need some preliminaries. Consider the problem of 
finding solutions of the form ip(r,uj) = R(r)fj P (u>) of equation with v = 
and D = B d {0, 1), where d > 2. We recall that: 

A = 7|-^ + ( d - 1 )^ 1 |- + r_2A ^- 1 ' 

(pry or 
where A S d-i is Laplace-Beltrami operator on S l , 

&s*-ifj P = -j(j + d-2)f jp . (5.2) 
Then we obtain the following equation for R(r): 

_ R » _ *Z±b> + j{j+ r d 2 - 2) R = Si?. (5.3) 
Taking R{r) — r ~ R(r), we get 

r 2 R" + rR! + \Er 2 - (j + ^—^-) ] R = 0. (5.4) 



This equation is known as the Bessel equation. For E = k 2 ^ it has two 
linearly independent solutions Jj + i=2.(kr) and Yj | <i-2 (for), where 

^— ' I (ra+ 1)1 m + a + 1 

m— 7 

J a (^) cos ira - J_ Q (2;) 

Ya(z) = : for a f Z, (5.6) 

sm 7ra 

and 

Y a (z) = lim Y a ,(z) for aeZ. (5.7) 

We recall also that the system of functions 

{ip jp (r,u) = Rj(k, r)f ip {u) : j € NU {0}, 1 < p < Pj } , 
is complete orthogonal system (in the sense of L 2 ) in the space 
of solutions of equation (TET]) in D' = B(0, 1) \ B{0, 1/3) 
with v = 0, E — k 2 and boundary condition ^| r= i = 0, 



(5.8) 



where 

Rj{k,r) = r-^r (Y j+ ^{kr)J j+i ^(k) - J j+ ^{kr)Y j+ ^ (fc)). (5.9) 
For the proof of (|5.8[) see, for example, |S]. 

Lemma 5.1. For any p > 0, integers d > 2, n > 10(p+ l) 2 and zeC, |z| < p, 
the following inequalities hold: 

1^/2^ 3^/2)^ 

2r^TT) - |Jq(z)I - 2rfo + i)' (5 - 10) 
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iJ4ooi<3 (| *jfflr \ ( 5 - u ) 

r(a) 

±-(\ z \/2)- a r(a) < \Y a (z)\ < A(|^|/2)-«r(a) (5.12) 
|yi(2r)|<-(M/2)— ^(a + l) (5.13) 

7T 

where ' denotes derivation with respect to z, a = n+^ 2 - and T(x) is the Gamma 
function. 

In fact, the proof of Lemma ISTTl is given in [5] (see Lemma 3.3 of [5]). It was 
shown in [8] that inequalities (|5.10[) - (|5.13[) hold for any n > no, where no is 
such that 

n > 3, 

exp ( ) - 1< 1/2, 

y \no + l) ~ 1 ' (5.14) 

max(l,(p/2) 2 »"+ 1 ) p 2 (p/2) 2 "°e^ 2 / 4 

3?r TV ^ h O 2 H TV \ - - 1 / 2 ' 

r(n ) 2n -p z T(n ) 

(see formula (6.18) of 0). The only thing to check is that n = [10(p+ l) 2 ] - 1 
satisfy (|5.14l) . where [•] denotes the integer part of a number, The first two 
inequalities are obvious. The third follows from the estimate 

1 \ no— 1 
n - 1 x 



r(no) - (n - 1)! > ) ■ ( 5 - 15 ) 

The final part of the proof of Lemma \'3 . 2 1 consist s of the following: first, we 
consider the case when E = k 2 ^ and 

ji = max{ 7 i,j 2 } > 10(1 + \k\) 2 . (5.16) 

Let ipi, ip2 denote the solutions of equation (11.11) with boundary condition 
ip\dD = fj 2 p2 an d potentials v% and t>2, respectively. Using Lemma 3.1 for 
V\ and V2, we get that 

Ui ~ fhWv^D) < 2(l + (N + \E\)Q) , (5.17) 

where 

Q = IK-A + ui - E)-% 2{D) ^2 {D) + \\{-A + v 2 -E)- 1 \\ h 2 {D) ^2 (D) , (5.18) 

Note that ipi -ip 2 is the solution of equation (JOJ in £>' = B(0, l)\B(0, 1/3) 
with potential w = and boundary condition V'|r=i = 0. According to (|5.8[) . we 
have that 

tpt -^2=^ c iv^iv m D ' (5.19) 



for some Cj p , where 



i )jp (r,Lj)=R j (k,r)f jp (u). (5.20) 
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Since Rj(k, 1) = 0, we find that 



dRj{k,r) 



dr 



d (r d i 2 Rj(k, r) 



r=l 



dr 



(5.21) 



For j>10(l + |fc|) 2 , using Lemma T5. 11 we have that 



dRi(k,r) 




dr 


r=l 



Y a (k)J a {k) 



\k\ 



Y a (k) J a {k) 



< 



f (\k\/2)-^T( a + l) (\k\/2)°-iT( a + l) 
\ (|*|/2)-<T(a) (|fc|/2)«r(a) 



(5.22) 



< 6|Jfe| 

lk"^ y a( fcr )llL2({l/3<| 2: |<2/5}) N 2 



6a, 



> 



\Y«(k)\ 

2/5 (i (\k\r/2)- a r( a y 1 



1/3 V3 (|fc|/2)-«r(a) 



> 



(5.23) 



\r ^ J a (fcr)|| )L 3({i/3<| a! | < 2/5}) 



< 



\J a (k)\ 
2 / 5 f 3 (\k\r/2rT(a) 



< 



1/3 V (|fc|/2) Q I» 



' (3(2/5)«) 2 r 



(5.24) 



where a = j + Note that if j > 10(1 + \k\) 2 then j + 4=Z > 3. Combining 
(j5T2"3"j) and (|5T2"3]> . we get that 

x|<2/5}) 



|Fa(fc)J Q (fc)| 



> 



> 



2 1\ 1 



5 3/3 



1/2 



-<•> ^" -3(2/r,r ) > i^q( 5 /2) q 



(5.25) 



dR 3 (k,r) 




dr 





Combining (pT2"2l and (|5T2"5]) . we find that 

< 1000«(5/2)-«||^p(^)lk 2( {i/3<|c C |<i»- (5.26) 
Proceeding from (|5.19[) and using the Cauchy-Schwarz inequality, we get that 



2 ({1/3<|z|<1}) 



IIV'jp(£ , )IIl2({i/3<| 2 ;|<i}) 

Using (|5.19[) . we find that 

[f jlP1 ,(&i(E) f hp ,) = 
dR n (k,r) 



< IMg) -^0(^)11^(5(0,1)) , g 2^ 



>(^)IIl2({1/3<|^|<1}) 



9(^1 - lp 2 ) 



fj 



dr 



fi 



dv 3D 

dRj 1 (k, r) 



dr 



(5.28) 



r=l 
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Combining lpB| , ([5T2"6"|) . (j5~2"7| and we obtain that 



(/ J1P1 , - * 2 (£)) / J2P2 ) < C(d)2^\\^ - ^11^(5(0,!))- (5.29) 

Combining ([5~T7|) and (pT2l?l) . we get (|3~T3)) for ji_> j 2 and E ^ 0. 

For ji < j 2 we use the fact that 4>*(2£) = &v(E) in order to swap ji and ji, 
where $* denotes the adjoint operator to Thus we complete the proof of 
Lemma 13.21 for the non-zero energy case. 

Estimate (|3 . 1 3|) for the zero energy case follows from Lemma 1 of [16] . 
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